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Abstract 

We study an integrable model of one-dimensional strongly correlated electrons at fi- 
nite temperature by explicit calculation of the correlation lengths of various correlation 
functions. The model is invariant with respect to the quantum superalgebra f7g(s[(2|l)) 
and characterized by the Hubbard interaction, correlated hopping and pair-hopping terms. 
Using the integrability, the graded quantum transfer matrix is constructed. From the an- 
alyticity of its eigenvalues, a closed set of non-linear integral equations is derived which 
describe the thermodynamical quantities and the finite temperature correlations. The re- 
sults show a crossover from a regime with dominating density-density correlations to a 
regime with dominating superconducting pair correlations. Analytical calculations in the 
low temperature limit are also discussed. 

PACS: 75.10.Jm, 05.50-fq 
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1 Introduction 

In recent decades, integrable models in one-dimensional strongly correlated electron systems 
have attracted considerable attention in relation to high- Tc superconductivity. In general, 
integrable models are solvable by means of the Bethe ansatz. The critical exponents of the 
correlation functions depend only on the i?-matrix of the underlying infinite dimensional sym- 
metry and the geometry of the Fermi sea of the ground state (see the book and references 
therein). This important fact shows that all integrable models with a given ii-matrix exhibit 
the same universal critical behavior. 

On the other hand (l-l-l)-dimensional critical phenomena are connected with conformal 
field theory (CFT) . The central charge and conformal dimensions are calculated by the scal- 
ing behavior of the ground state and the low-lying excitations The conformal dimensions 
determine the long distance asymptotics of all correlation functions of any local operators. For 
integrable models, one obtains the conformal dimensions by use of the Bethe ansatz in terms 
of the so called dressed charge function satisfying a linear integral equation. Thus the Bethe 
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ansatz together with CFT provide a powerful tool to study the correlation functions especially 
at the critical point T = |Q,|5|. 

For finite temperatures however, one encounters technical difficulties. Despite the fact that 
the models are still solvable even off criticality, the above powerful methods lose effect due to 
the collapse of conformal invariance. Qualitatively, we see the long-distance behavior of the 
correlation functions show exponential decay and calculate the low temperature asymptotics by 
extension of the conformal mapping (see for example Q ) . However we need another approach 
for the quantitative understanding of the correlation functions at finite temperatures. 

As a totally different approach, the quantum transfer matrix (QTM) has been proposed 
recently to overcome such difficulties [^-0]. Originally the QTM method has been developed 
as an alternative method to the traditional thermodynamical studies |p!7|-|T9|. Utilizing the 
Trotter mapping, one deals with a two-dimensional classical system instead of the original 
one-dimensional quantum system. Using the underlying integrability, one constructs a family 
of commuting QTMs ]T5| , . The commuting QTMs are diagonalized by means of the Bethe 
ansatz. Thermodynamical quantities are expressed in terms of the sole largest eigenvalue. 
Also, this procedure has remarkable advantages for the studies of correlation functions, which 
is the main topic of our paper. Explicitly, the correlation lengths can be calculated by taking 
ratios of the largest eigenvalue and the sub-leading ones 20-23|. 



In this paper we study a strongly correlated electron system using the approach sketched 
above. As a result, we evaluate the six crucial correlation lengths for arbitrary particle densities 
and finite temperatures: (i) the longitudinal and (ii) transversal spin-spin correlations, (iii) the 
density-density correlations, (iv) the one-particle Green function, (v) the singlet and (vi) triplet 
superconducting pair correlations. These are the ffist quantitative calculations for correlation 
functions of strongly correlated electron systems off criticality. 

The Hamiltonian of the model |24| on a periodic lattice of size L is defined by 



L 



where cj^ (cj,o-) denotes the fermion creation (annihilation) operator at jth site satisfying the 
canonical aniticommutation relations 

{4,^' 4,t} = {Cj>> Cfc,r} = 0, {c] Ck^r) = SjkSar, (1-2) 

and Uj^cr is the number operator, i.e. nj^cr = cjg.Cj,o-. The ground state properties and the 
long-distance behavior of the correlation functions have been studied in |24| by the Bethe 
ansatz and CFT. Due to the massive spin excitations (7 7^ 0), the ground state properties 
are described by a one-component Tomonaga-Luttinger liquid (Luther-Emery liquid). For 
this system the one-particle Green function decays exponentially fast, the density-density and 
certain pair correlations are quasi-longranged with algebraic decay. As the model is one- 
dimensional it does not exhibit off-diagonal long-range order. However, the crossover behavior 
from dominant density-density correlations to dominant (singlet) pair correlations driven by 
the density of particles manifests superconducting properties. More quantitatively, on the basis 
of scaling arguments and RG concepts we would expect that higher dimensional interactions 
such as interchain couplings are relevant and may lead to true condensation of Cooper pairs. 

The model is related to a trigonometric i?-matrix which is a solution to the graded YBE 
associated with the four-dimensional irreducible representation of the quantum superalgebra 
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C/y(s((2|l)) (see also |27| for the rational case 7 = 0). The diagonalization of the Hamilto- 
nian for general 7] and 7 was done in and the diagonalization of the associated row-to-row 
transfer matrix was done for the rational case 7 = in [28, 2^]. The thermodynamics have been 



investigated by the traditional thermodynamic Bethe ansatz based on the string hypothesis 



in 1 3C ] and the present approach in . 

The asymptotic behavior of correlation functions of this model have been studied so far only 
for the ground state where a crossover was established from a normal to a "superconducting" 
regime driven by the particle density. Here we want to extend this study to finite temper- 
atures where we are especially interested in the extension of the "superconducting regime" 
in the density-temperature phase diagram. For this purpose we define the QTM acting on a 
Z2-graded vector space. As a consequence, the genuine fermionic statistics are built into the 



QTM algebraically through the Grassmann parity p[j] (see eq. ( p. 26 )). Regarding the ther- 



modynamical quantities, one might content oneself with the corresponding null Grassmann 
parity modelQ because the fermionic statistics does not affect the thermodynamical quantities 
directly. However, one must be careful to keep track of the fermionic statistics when consid- 
ering the physical quantities changing the particle number. As the auxiliary and quantum 
space are the same graded vector space, we easily obtain the eigenvalues of the QTM by using 
the results from the algebraic Bethe ansatz for the row-to-row transfer matrix. The resultant 
eigenvalues can be expressed in the dressed vacuum form (DVF) via the Bethe ansatz equation 
(BAE). The difference between the DVF of the QTM and the one of the row-to-row transfer 
matrix consists only of their vacuum parts. Hence, in the actual calculation to the DVF, one 
utilizes the methods developed in the finite-size correction problem |3^. In consequence, all 
information of the BAE is contained in three non-linear integral equations (NLIEs), which are 
valid for a system with any Trotter number A^. In this approach, the limit N ^ 00 can be 
taken analytically. This finite set of NLIEs allows for highly accurate numerical calculations of 
physical quantities such as the above mentioned correlation lengths. In particular, we observe 
temperature dependent crossover phenomena between regimes dominated by density-density 
or (singlet) superconducting pair correlations, respectively. 

At low temperatures, the three NLIEs are reduced to only one linear integral equation, 
which reflects the one-component Tomonaga-Luttinger (Luther-Emery) liquid properties. This 
linear equation is evaluated analytically and reproduces the predictions from GET. The for- 
mulations in this paper are applicable to various kinds of strongly correlated electron systems 
which have superalgebra invariance, for example the supersymmetric extended Hubbard model 



(s[(2|2)-invariance) |3|,p4|. 



The layout of this paper is as follows. In section 2 we define the commuting QTM con- 
structed by the graded i?-matrix. The QTM is diagonalized by the algebraic Bethe ansatz in 
analogy to the ordinary row-to-row transfer matrix. In section 3 we study the eigenvalues by 
introducing certain auxiliary functions. Investigating their analyticity, we determine the NLIE 
for the largest and sub- leading eigenvalues. We obtain six sets of NLIEs characterizing the 
above mentioned correlations. The numerical results for these NLIE are discussed. Section 4 
is devoted to the analytical calculation of the low temperature asymptotics of the correlations 
characterized by gapless excitations. Section 5 contains our conclusions. Details of technical 
calculations are summarized in various appendices. 



^ We obtain the null-parity model satisfying the ordinary Yang-Baxter equation by multiplying the original 
graded i?-matrix by suitable minus signs. 
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2 The graded quantum transfer matrix 



The model (^]^) is solvable under the conditions 

Here the physical region of the two parameters a and 7 is restricted to 7 > and q > 
(repulsive) or a < — 1 (attractive). [This does not include the pure Hubbard model which 
corresponds to 7/ = 7 = ip = and U ^ violating ( |2.1| ). Note however, that a unified 
treatment of the C/q(s[(2|l))-invariant generalized Hubbard model and the standard Hubbard 



model might be possible on the basis of the work in |35].] 

Several prominent models are comprised as special cases: the super symmetric t — J model 
(7 = 0, a — > 0), the pure correlated-hopping model (|a| —>■ 00) and a model of hard core 
particles of three species (doubly occupied site, single electron with spin pointing either up or 
down) (q — !)• For any finite 7 the spinon excitations have a mass gap and so do single 
particle excitations. The pair excitations (singlet) and the particle-hole type excitations are 
massless. The former excitations become gapless for the repulsive case (a > 0) in the rational 
limit 7 = 0. Consequently for 7 > the one-particle Green function shows exponential decay 
and hence the momentum distribution function has no singularity at the Fermi point. 

The classical counterpart of this model is given by a trigonometric i?-matrix which is 
associated with the one parameter family of the 4-dimensional irreducible representation of 
the quantum superalgebra Uq{sl{2\l)) [26| (see also ref. for the rational case 7 = 0). 
Choosing the following basis 

|l)^|n) = c|c{|0), |2)^U) = cI|0), |3)^|T) = c{|0), |4) = |0), (2.2) 

with grading |4) even (bosonic) and |2), |3) odd (fermionic), one determines the following 
36 non-zero matrix elements as 

n.,x[2a^ 22. s _ „33. ,N _ 1 44 _ [2a + 2 + 1.] 

>12/„,\ _ d13/„,\ _ d21/„,a _ d31 



Riiiv) = RiUv) = mHv) = Rtiiv) 



[2a + v] ' 

Rltiv) = Ritiv) = RtUv) = Rtliv) 



[2a + 2-vy 



[2a + v][2a + 2 - v] ' ' [2a + v][2a + 2 



[2a + v] ' ' ' [2a + vy 

41 e-^- [2a] [2a + 2] ,4 e^- [2a] [2a + 2] 

^14 J = — : — TT^ — -T, T' Rai\^) ~ 



[2a + v][2a + 2 - v] ' [2a + v][2a + 2 - v] 



.32. ^ _ -e^"/^[2][^]-[2a][2a + 2] .3 _ -e-^-/^[2][v] - [2a][2a + 2] 
[2a + v][2a + 2-v] ' ^32 (^'J - [2a + v][2a + 2 - v] ^ 
42, ^ _ p43. _ e-^-/2[2a + 2] .4 , , _ .4 , , _ e^^/^[2a + 2] 



Rlliv) = -Riliv 



[2a + 2-t;] ' ' ' [2a^2-v\ 

^[2a][2a + 2]e-^(i+^)/2 [v\ 



[2a -h v\ [2a + 2 



Rlliv) = -Rtliv) 



V[2a][2a + 2]e^(^-^)/2[7;] 
[2a + v][2a + 2 - v] 
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>l4,.,^_ z,32,.,^_ V[2«][2« + 2]e^(i+'')/2H 



[2a + v][2a + 2- v] 



. _flg(„) . (2.3) 



[2a + v][2a + 2 - v] 

a/3 

,75 



where -^^^(w) are defined by 



i?(^)[a)®|/3) = |7)^|<^>i2r/.(«) (2.4) 



with [v] := sh(7u/2) and the graded tensor product satisfies P\a) (g)^ := (— )pNp[^1|/J) (g)^ [a) 
where P denotes permutation. This i?-matrix is known to satisfy the graded Yang-Baxter 
equation (GYBE) pl 



i?i2(n - v)Ri3{u)R23iv) = R23{v)Ri3{u)Ri2{u - v), (2.5) 
which exphcitly reads in terms of the matrix elements 

_ ^)^71A(^)^7273^(^)(_){P["1]+P[71])P[72] 

= R^4^l{v)R^lZliu)Rl\l\{u - t;)(-)(P[/3il+P[7i])P[72]^ (2.6) 

where the Grassmann parity p[j] is defined by p[l] = p[4] = (bosonic) p[2] = p[3] = 1 
(fermionic). In eq. (|2.6| ) and later we sum over repeated indices. Note that the parity of the 
i2-matrix is even, i.e. 

p[R] = p[a] + p[P] + p[-f] + p[6] = mod 2. (2.7) 
To relate the model ( [1.1| ) to the above classical model, we construct the graded transfer matrix 

T{v) = StlalRaLiv) . . . Ra2iv)Raliv)} . (2.8) 

Using the i?-matrix parity (^]3), one writes explicitly the transfer matrix by means of the 
matrix elements ( ^.3]) 

T^lt.iii^) = (-)PW+S'-(P["^l+^[^^])S^^-^[^^-lii^,_^f^(^) . . . R;f^liv)R:t{v). (2.9) 



Due to the GYBE ( |2.5D , the transfer matrices T{v) commute for different spectral parameters 
V and v' 

[T{v),T{v')]=0. (2.10) 
Using the fact that the local Hamiltonian Ti-jj+i can be expressed as the derivative of Rjj j^i{v) 

2sh(a + 1)7 „ d 

J^jj+i 7 ^jj+i-jz^jj+n^) 

+ 2ch(a + 1)7 - 2(nj_| + Uj^i + nj+i^-f + nj+i_;)ch(a + 1)7, (2-11) 
we can expand the transfer matrix T{v) with respect to v as 
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Here we introduced the superpermutation operator Pjj+i whose matrix elements are 

[PmiVi = [Rn+mVjp = (2-13) 

Note that H' denotes the Hamiltonian (LI) with shifted chemical potential and ground state 
energy 

n' = n- 2Lch(Q + 1)7 + 2A4ch(a + 1)7, (2.14) 

where Afe is the electron number operator. 

In order to treat the finite temperature case, we consider another i?-matrix defined by 

Wjf^iv) = {-)Pif^M'']+Ph])Rp^{v). (2.15) 

This matrix R satisfies the following type of GYBE 

Ruiv - u)Rni-u)R23iv) = R23iv)Ri3{-u)Ri2iv - u), (2.16) 

which reads explicitly 

= <%^|(^^)<l:(-n)i2l7i,(. - n)(-)(^[/^^l+^[^^])^[^< (2.17) 
Using this, one constructs the transfer matrix 

T{v) = StTa{RaL{v) . . .Ra2iv)Raliv)} , (2.18) 

explicitly reading 

TZti>) = (-)-W+St^^(p[".]+P[ftl)EUp[/^.]7^;_^f^(,) . ..^fj^iv)R:f^{v). (2.19) 
We expand this transfer matrix in the same manner as ( 2.12| ), 

Combining the above relations and using T(0)r(0) = 1 (note that T(0) and r(0) denote the 
right and left shift operators, respectively), one obtains 

T{n)T{n) = 1 - H' + 0{v'). (2.21) 

sh(a + 1)7 

From the above equation, it follows that 

eM-(3n') = hm [T{un)TM]''/\ ^ 2/3sh(« + 1)7 ^ 

Here (3 is the reciprocal temperature; /3 = 1/T and the height N of the fictitious underlying 
square lattice is even and referred to as the Trotter number. Hence the free energy / of the 



original quantum system (LI) is given by 

lim lim —— 

L— ►00 N^oo Ljj 



f = - lim lim -i- In (Tr [T{un)T{un)] ^^^) + 2(1 - ne)ch(a + 1)7, (2.23) 
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where Ue denotes the particle density of the system. Calculating Tt[T{uN)T{uN)]^^'^ from the 
eigenvalues of T{u]\f)T{u]\f) is a serious problem, because the spectrum is infinitely degenerate 
in the limit of infinite Trotter number un^O (for oo). To avoid this difficulty, we 

transform Tr[T(njv)T(uAr)]^/^ as 

N/2 

Tr Yl StIa2^,a2k-l{Ra2kL{uN) • • • Ra2kl{uN)Ra2k-iL{uN) ■ ■ ■ Ra2k-il{uN)} 
k=l 

L N/2 

= StTYlTTjl[Ra,,,{uN)Ra2,^,j{uN). (2.24) 
j=l k=l 

Now we introduce the QTM, 

N/2 

Tq{v) = Trj Ra^k^N + v)Ra2^_^j{uN - v), (2.25) 
k=l 

explicitly 

J. /3i/32.../3iv ^ (_)Ef=2{pN]+P[/3j])Eillp[A] 

X R^f-y.-. i^N + v)Rt--\Z--. (nN-v)... RiZ i^N + v)Ri\l- {um - v). 

(2.26) 

Due to the GYBE (|2j) and ( |2T6D , the QTM is commutative, i.e. 

[Tq{un,v),Tq{un,v')\=Q. (2.27) 

We will see that the largest eigenvalue is separated by a gap from the rest of the spectrum for 
any persisting in the limit ^ oo as long as T > 0. Since the two limits in eq. ( p.23| ) can 
be interchanged, for a proof cf. [^,|7|, we take the limit L ^ oo first. As there exists a finite 
gap between the largest eigenvalue A™^^(0) and the sub- leading eigenvalues A*''^'^(0), we can 
write 

/ = - lim ^lnA'^^^(0) + 2(nc-l)ch(a + l)7. (2.28) 

N^oo p 

In this approach, the thermodynamical completeness — lim^^o = In 4 follows easily from 
limiv-»oo A™'^^(0) = 4 which is obvious from the definition of R{v) ( |2l^ ) and R{v) ( |2.15|) . 

Most significantly the present method makes it possible to calculate various correlation 
lengths at finite temperature through the relation 



— = — lim In 

^ Af-+cx> 



^max(o) 



(2.29) 



The QTM can be diagonalized by means of the algebraic Bethe ansatz method. Due to the 
fact that auxiliary and quantum space are the same Z2-graded vector space, one can utilize 
the results for the algebraic Bethe ansatz for the row-to-row transfer matrix case p8|, ^p|. 
Replacing the parameters un and v by 

7 . 7 /?sh(a 4- 1)7 „„^ 
-V tv, -UN ^UN = j;^ , (2.30) 



^ We remark that the formulation in is based on the non-graded rational (7 = 0) 7?-matrix which 
is connected to our graded case by _R^^ — » ( — ^Phl^l'^li^^* . The basic difference between the standard and 
graded formulation lies in the existence of extra phase factors in the dressed vacuum form and the Bethe ansatz 
equation. 
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and choosing the state (|1) |4))®= as the pseudo vacuum state (reference state), one 
writes the eigenvalue A{v) of the QTM in the dressed vacuum form (DVF) 

A(.) = ^^(,)^i(^ + §7(2« + l))^,,, ^ ^^^^^g,(. + |7(2a + l))..(. + n)^,, 
Qi{v + ^l) gi(w + §7) q2{v) 

+ Mv) j— r-. — e^"' + (p3[v) — , (2.31) 

where the functions (j)i{v), (j)2{v) and (psiv) are defined by 



(j)i{v) 



N 



sh{iv + UN)sh{iv + un — ^)s\i{iv — un + 07) ^ ^ 

sh{iv + ujy — (a + l)7)sh(it; + un + a7)sh(iT; — ujy — 07) 



JV 



\sh(«f — — Q;7jsh(«f + UN — (a + 1)7) / 

= 0i(-t')U^-Q-i, 

and the chemical potential /u (shifted hy fi + 2ch(a + 1)7 fi, see Appendix A) have been 
introduced. The functions qi{v) and (72 (w) are written in the form 

n m 

9i(«) = ~ "^i^^^' 92 (^^) = sin(i; - vf^), (2.33) 

where the unknown parameters {wj"'^^}je{i,2,- - ,n} {^j^^}jG{i,2, -- ,m} ^-^'s the Bethe ansatz 
rapidities determined from the Bethe ansatz equation (BAE) 

92(wf ^ - It) Mvf^ - It) ' qi{vf - |t) 92(?^f ^ - n) 

At the end of this section, we emphasize that our QTM is based on the graded formulations, 
which means the genuine fermionic statistics of the model are properly built into the algebraic 
structure of the QTM through the Grassmann parity p[j] (see ( p. 26 )). Nevertheless the DVF 
and the BAE have quite simple forms. This indicates that the differences between the null and 
non-null Grassmann parity models in one dimension are embedded into differences of boundary 
conditions. 



3 Non-linear integral equations 



In this section we analyze the eigenvalue A{v) in DVF form (2.31) by selecting certain auxiliary 
functions including the information of the BAE. From now on we consider the repulsive case 
a > 0. For the largest eigenvalue, the auxiliary functions and NLIE have been introduced 



in |31|. For some special eigenvalues with certain symmetries including the largest eigenvalue, 
one needs only two auxiliary functions |3l|| . 

In order to treat the general case, we introduce the following three auxiliary functions. 



Ai(x)(A3(2;) + A4(x)) 



A2(2;)(Ai(x) + A2(2;) + A3(x) + A4(2;)) 

A2(x) 

4 fx) 



Xsix) + A4(x) 
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ai{v) 



Xi{x 



2ti(v) = l + ai{v). 



(3.1) 



jQ7 



where the functions Xj{x) denote the jth term on the right hand side in the DVF (2.31). 

These auxihary functions satisfy certain functional relations. Exploring their analyticity, 
one transforms these functional relations into a closed set of NLIE. Here we consider the NLIE 
for the largest eigenvalue and some sub-leading ones which describe the correlations. The 
details of the calculations are deferred to Appendices B and C. 



3.1 Largest eigenvalues 

First we consider the largest eigenvalue. As described above, one finds the NLIE for this 
eigenvalue in |31]. For completeness we study this case by using the above three auxiliary 
functions. 

The largest eigenvalue belongs to the sector n = N and m = N/2 in the BAE ( |2.34 ). 
From numerical calculations with finite Trotter number A^, we observe the above auxiliary 
functions have the following analytical properties: The function 2to(f) (2lo(u)) is analytic and 
non-zero in a finite strip in the lower (upper) half plane —7/2 < \mv < (0 < Imu < 7/2)^, 
and the function 2li(f) is analytic and non-zero along the real axis. Hereafter we call the 
region —7/2 < Imu < 7/2 "physical strip". Due to this analyticity, we can transform the 
functional relations satisfied by the above auxiliary functions into NLIEs by using the Fourier 
transform together with Cauchy's theorem. In this procedure one takes the Trotter limit 
N ^ CO analytically (see Appendix B). Consequently, we obtain the following closed set of 
NLIE. 



Inao(v) = I3il){v) + k* ln2to(i') + k* \n%i{v) + 
Inao(w) = I3i){-v) + A: * ln2to(t;) + fc*ln2li(u) + Pn, 

lnai(?;) = I3{ij{v) + iIj{-v)) + A:*ln2lo(u) + k*\u.%o{v) + ki *ln2ti(i;) + 2/?//, 



(3.2) 



with the kernels and leading terms 

sh7 



k{v) 



2s\i{iv)s\i{iv — 7) ' 

sh^(a + 1)7 



k{v) = k{—v), ki{v) = k{v) + k{v), 



sli{iv + ^aj)sh{iv — ^7(0 + 2)) ' 
where the symbol * denotes the convolution defined by 



/ * g{v) = - [ f{v- x)g{x)dx = - [ f{x)g{ 

Jc ^ Jv-C 



V — x)dx. 



(3.3) 



(3.4) 



The integration contours C for convolutions with InSlo, InSlo and ln2li in ( p.2| ) should be 
taken by a straight line ranging over a full period vr with imaginary part —6, +6 and {5 is 
arbitrary but fixed in the range < (5 < 7/2), respectively. Through the solution to the above 
NLIE, the largest eigenvalue is expressed by 



In A" 



[V] 



■^{v) + C * ln2lo(?;) + C * ln2to(w) + (C + C) * ln2ti(u) 



(3.5) 



In fact these auxiliary functions have trivial zeros and poles of order N/2 which come from the vacuum 
parts in the DVF. These zeros and poles partly appear in the above "analyticity strips" and result into the 
leading terms of the NLIEs. 
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where the kernels C(^)) C(^) the leading term ^(v) are defined by 



sh(zf + (q + l)'y)sh[iv — [a + 1)7) 2sh(a + 1)7 



From eq. (|3.2D , we find ao{v) = ao{—v) and the real function ai{v) is symmetric with respect 
to both real and imaginary axis. Hence the threefold set of NLIEs ( |3.2| ) can be reduced to 



only two NLIEs which are identical to the ones in [31 



3.2 General NLIE for arbitrary eigenvalues 

The sub-leading eigenvalues of the QTM ( p.31| ) are treated in a way similar to the largest eigen- 
value. The main difference lies in the existence of additional zeros and poles of the auxiliary 
functions. Therefore one observes additional terms in the NLIE ("sin" -terms). Consequently 
all correlation functions are characterized by the distribution patterns of these additional zeros 
or poles. 

Often in the case of excited states the second factors in the convolutions are logarithms of 
periodic functions, however with non-zero winding numbers, ln2l(7r/2) — ln2l(— 7r/2) = 2mTi 
[n ^li). This gives rise to a second set of additional terms in the NLIE ("cos" -terms). 

Thus we obtain the following general NLIE with additional terms In 939 (^)) lii</?i(t') and 
In^oM- 

lnao(t;) = [iil){v) + k* ln2to(f) + k* ln2li(u) + hnpQ{v) + 

\nao{v) = ;3^(-w) + fe*ln2lo(v) + A?*ln2ti(u) +ln^o(u) (3.7) 
lnai(i;) = I3{ip{v) + V(-^)) + k* ln2lo(i') + k* hi%o{v) + ki* ln2li(f) + ln(/7i(v) + 2/3/x, 

The additional terms in the above general NLIE must be determined for each eigenvalue 
separately. Note that the relation (pi{v) = ipo{v)TpQ{v) is always valid due to certain algebraic 



structures of the auxiliary functions (see Appendices B and C). In (3.7) and later, we take 
the same integration contours as in the case of the largest eigenvalue. Thus the corresponding 
NLIE for arbitrary eigenvalues (including the largest eigenvalue) are written in the following 
general form: 

lnA(?;) = ^'(t;) + C*ln2lo(f) + C*ln2to(t;) + (C + C) * ln2li(w) + In x(w), (3.8) 

where the term \\ix{v) is also determined from the distribution pattern of the additional zeros 
and poles. 

As concrete examples, we determine the additional terms for the eigenvalues corresponding 
to the asymptotics of the: (i) one-particle Green function, (ii) transversal spin-spin correlations, 
(iii) density-density correlations or longitudinal spin-spin correlations, (iv) sub-dominant corre- 
lations of case (iii), (v) singlet and (vi) triplet superconducting correlations (see Appendix C). 



3.3 Numerical analysis of the NLIE 

Here we evaluate the above NLIE numerically. As seen in the above section, the NLIEs are 
exact expressions for the thermodynamical quantities and various correlations at any finite 
temperature. In contrast to the standard TEA they close at a finite level. Hence one obtains 
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physical quantities explicitly by a highly accurate numerical analysis. To keep the electron 
density constant, we adopt temperature dependent chemical potentials ^(T) determined from 



d(ne(T,/i(r))) 



dT 



d 

dT 



0. 



(3.9) 



In Fig. H we depict /u(T) for a = 1, 7 = 1 and various particle densities. We solve the NLIE 
by the iteration method. In each step, convolution parts are calculated by making use of the 
Fast Fourier transform. The subsidiary conditions which determine the zeros ^i) etc. 
for example, eq. ( C.5D ) are solved by the Newton method. 



see, 




Figure 1: Temperature dependence of the chemical potentials for interaction parameters a = 1, 
7 = 1 and various particle densities. The symbol * denotes the zero temperature chemical 



potentials derived by ( A.5 ) and ( A. 9 ) in Appendix A. In our units the Boltzmann constant is 
kB = 1. 



Let us consider the NLIE characterizing the correlation lengths via formula p^ ). Fig. 



shows the temperature dependent correlation lengths (multiplied by T) of the one-particle 
Green functions for various particle densities. The case a = 7 = 1 is used as a concrete 



example throughout this paper and is equivalent toU /2 = tp = 1 in (LI). Due to the finite gap 
of single particle excitations, we observe that the correlation length multiplied by temperature 
shows monotonous decay with decreasing temperature with zero limit for zero temperature. 
We observe a crossover behavior of the correlation lengths at low temperature close to particle 
density rig = 1.5: for < 1.5 (ne > 1.5) the values of ^ increase (decrease) with n. A similar 
behavior is observed in other correlations described by massive excitations. The transversal 
spin-spin correlations, the triplet superconducting pair correlations are depicted in Fig. ^a and 
in Fig. ^3, respectively. We also depict the temperature dependent oscillatory terms of the one- 
particle Green functions, transversal spin-spin correlations, and density-density correlations 
in Fig. |2|b. Fig. ^3 and Fig. |5|b, respectively. Regarding the one-particle Green functions, 
one observes clearly a temperature dependence of the "Fermi momentum" kp. In the low 
temperature limit T — > 0, the wave vector converges to the expected value kp = 7rne/2 which 
indicates the significance of the Fermi surface for one-particle excitations in the Tomonaga- 
Luttinger liquid at T = 0. As seen in Figjl] (and also expected) the chemical potential grows to 
-|-oo (—00) for T ^ 00 if the particle density is kept fixed at a high (low) value. This indicates 
a broadening of the momentum distribution at finite and in particular at high temperatures. 
Hence, processes of particle-hole type closer to the middle of the band become statistically 
dominant. Consequently, the wave vector kp decreases (increases) with increasing T if the 
T = value was large (small). As shown in Fig. ^b and Fig. ^3, oscillatory terms of the 
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Figure 2: (a): Temperature dependence of the scaled correlation lengths T ■ ^(T) for the one- 
particle Green function, (b): Temperature dependence of the wave vector of the oscillation 
terms. At zero temperature one observes the ^-oscillation (kp = 7rne/2 depicted by the 
symbol *) characteristic for the one-particle Green function. 
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Figure 3: (a): Temperature dependence of the correlation lengths of the transversal spin-spin 
correlations, (b): The wave vector 2A;f of the oscillations which is characteristic for the spin- 
spin correlations. At zero temperature one observes the strict relation of wave vector with 
particle density 2kp = irue (analytical values are depicted by the symbol *). 
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Figure 4: Temperature dependence of the correlation lengths for singlet superconducting pair 
correlations (a) and triplet pair correlations (b). The symbol * in (a) denotes results of 
the analytical calculation in the low temperature limit presented in section 4 confirming the 
predictions from CFT in Appendix A. 

density-density and transversal spin-spin correlations converge to 2kp vrne- This behavior 
corresponds to the excitations carrying a 2kF momentum. For the same reason as mentioned 
in the case of the one-particle Green function, this 2A;f is shifted with increasing temperature. 
The quantitative difference of the oscillation terms for the transversal spin-spin correlations and 
the density-density correlations lies in the different nature of these excitations. The density- 
density correlations are described by particle-hole excitations from the left to the right Fermi 
point without changing the electron and spin numbers. On the other hand, the transversal 
spin-spin correlations are characterized by massive excitations changing the spin number. 

Fig. ^ and Fig. ^ show the density-density and the singlet superconducting pair corre- 
lation lengths, respectively. Due to massless excitations, the scaled correlation lengths T^{T) 
are finite at zero temperature. The low temperature asymptotics agree with the analytic 
calculations in section 4 and the prediction from CFT in Appendix A (see Fig. p4| ). At low 
temperatures one observes a crossover from dominant density-density correlations to dominant 
singlet superconducting pair correlations driven by the particle density (see also Fig. ^ in Ap- 
pendix A): below (above) a certain critical density pc, the dynamics is dominated by singlet 
superconducting pair correlations (density-density correlations). This behavior is changed at 
higher temperatures. Here at low (high) density p the density-density correlations (singlet 
superconducting pair correlations) dominate, see Fig. ^. Of course, these findings for high par- 
ticle density and temperature do not imply any tendency towards superconductivity. It should 
be born in mind that the overall length scale of the density-density and singlet superconducting 
pair correlations is rather short in comparison to the lattice constant. Furthermore, neither 
of these correlations dominate the dynamics with increasing temperature as the one-particle 
Green function takes over. This implies a picture of rather uncorrelated one-particle dynamics 
at higher temperature. 

Finally we comment on the sub-dominant density-density (longitudinal spin-spin) correla- 
tions in Fig. ^. The correlations are characterized by simple particle-hole excitations at each 
Fermi point. Correspondingly, the correlation lengths (multiplied by temperature T) approach 
the Fermi velocities fF/(27r) in the low temperature limit. 
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ii.=1.50 - 



Figure 5: (a): Temperature dependence of the correlation lengths for the density-density 
correlations. The wave vector of the oscillation term is depicted in (b). The symbol * denotes 
results of the analytical calculation in section 4 in the low temperature limit and the predictions 
from CFT in Appendix A. 

4 Low-temperature asymptotics 

In this section we consider the model at low temperatures. The largest eigenvalue and those 
sub-leading ones described by massless excitations are calculated approximately up to order 
0(1/ P). As a consequence, the NLIEs are reduced to only one linear integral equation and 
this linear integral equation is connected directly with the dressed functions (see Appendix A). 

Due to the massive spin excitations, the two auxiliary functions In ao{v) and In ao(f ) can be 
neglected as ao,ao e~^^ (e > 0). Hence the three NLIEs characterizing the grand potential 
and correlation functions reduce to the single integral equation 

Inai(u) = -f3e^°\v) + ki ln2li(t;) + ln(^i(i;), (4.1) 

where 

e(o){v) = -tpiv) - ip{-v) - 2fi. (4.2) 

and (fii is a collection of only the "cos"-terms of the function ipi listed in Appendix C. The 
"sin" -terms do not show up explicitly as they are taken care of by the modified contour Ci (see 
Appendix C). Throughout this section the convolution / * g{v) is understood in the manner 

f*9{v)-=- I f{v-x)g{x)dx. (4.3) 
71" Jc 

Numerically, one observes that the function ai{v) has a crossover behavior from ai{v) ^ 1 to 
Oi(f) ^ 1 (see Fig. |8|). The Fermi point Ap can be defined by the following conditions satisfied 
by Inai(v) for the largest eigenvalue (Ci is taken by a straight line on the real axis), 

lnai(±AF) = (1//3<1), (4.4) 

and for the excited states we demand Re Inai(ibAF) = for itAp on the real axis. 

4.1 0{P) and 0(1) approximation 

First we calculate order 0{P) and 0(1) approximations for the largest and those sub-leading 
ones characterized by gapless excitations: (i) particle-hole excitations at each of the left and 
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Figure 6: (a) Low particle density: crossover from dominant (superconducting) singlet pair 
correlations at low temperature to dominant density-density correlations at high temperature, 
(b) High particle density: opposite to scenario in (a) for density n = 1.25 (upper two curves). 
For n = 1.0 the crossing lies outside the shown temperature window. 




Figure 7: Temperature dependence of the sub-dominant density-density (longitudinal spin- 
spin) correlation lengths. At low temperature there is convergence to the Fermi velocity 
fF/(27r) (denoted by the symbol *). 

right Fermi points (define n+ and n~ as corresponding quantum numbers), (ii) particle- hole 
excitations of dc charges from the left to right Fermi point, but no change of the total charge 
number and total spin number, and (iii) singlet pair excitations near the Fermi points; an- 
nihilate or create ANc = 2nc {ric € Z) charges in a symmetric way near the Fermi points. 
The longitudinal spin-spin correlations {(Tjaf) and the density-density correlations {ujUi) are 
characterized by excitations (i) and (ii). On the other hand, the singlet superconducting pair 
correlations (cj+i^fCj^iCj+i^i-Cj^i) is described by (iii). 

(1) The largest eigenvalue 

In the low temperature limit, the function ln2li(f) can be written in terms of ai{v) as 
ln2li(i;) ~ Inai(f) for \v\ > Ap and In 211(11) ~ for \v\ < Ap (Fig. P). Hence we have 

In ai (v) = -/?e(°) {v) + h ^* In oi (v) + 0(1/0) , (4.5) 



15 




Figure 8: Crossover behavior of the functions Inai(f) and ln'Qli{v) (real parts) for = 7 = !, 
ne = 1.00 and P = 50. The Fermi point Ap defined by is found to be Ap = ±0.748. 



where the symbol * for the case of the largest eigenvalue denotes 



Ci 1 /" 

ki * Inai(f) = — / ki{x)lnai{v — x)dx 

= ~ij + J^j ki{x) In ai{v - x)dx. (4.6) 



Using the definition of the dressed energy e{v) ( |A.7| ), we find a solution to the above equation 

\nai{v) = -Peiv)+Oil/P). (4.7) 
Then we have an approximate value for Inai(f) at the Fermi points 

In ai(±AF) = + 0(1//?). (4.8) 
(2) Particle-hole excitations at each of left and right Fermi points 



This excitation is characterized by excited charges at the left (right) Fermi points — Ap (Ap). 
In the simplest, still characteristic case, one zero is separated by n~ {n~^) holes from the Fermi 
sea of zeros at the Fermi point =FAp. The holes (zeros) are to be circumvented by the deformed 
contour in anti-clockwise (clockwise) manner (see Appendix C). The function ln2li(T;) can be 
approximated in terms of ai{v) in the same manner as case (1). Hence we obtain the same 
approximate value as in ([4. 81) . 

(3) Particle-hole excitations from left to right Fermi points 

These excitations are described by dc charges moved from the left (right) Fermi point to 
right (left) one. This excitation entails a change of momentum 2kpdc. Correspondingly, the 
auxiliary function ^i{v) has 2dc zeros near the Fermi points with behavior different from those 
in the largest eigenvalue case: dc zeros (of "particle type") are distributed in the upper half 
plane near the left Fermi point and the remaining dc zeros (of "hole type") are located in the 
lower half plane near the right Fermi point, and vice versa. 

These zeros should be circumvented by the deformed contours in clockwise manner at the 
left Fermi point and anti-clockwise manner at the right Fermi point (see Appendix C). At low 
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Figure 9: Numerical results for the functions Im In ai{v) (a) and Im ln2ti(v) (b) for = 7 = !, 
Ue = 1.00, P = 50 and dc = 1. Due to the additional zeros (circumvented by the modified 
integration contour Ci), ln2li(?;) is discontinuous near the Fermi points Ap = ±0.748. 

temperatures, ln2li(w) is replaced by ln^i{v) ~ for \v\ < Af and ln2li(f) ~ lnoi(z;) + 27ridc 
for \v\ > Ap (see Fig. ^ for dc = 1). Thus we have the linear integral equation characterizing 
the excitation. 

Inoi(w) = -pe^°\v) + ki ln{ai{v) + 2TTidc) + 0{1/ /3), (4.9) 

where Ci denotes the integration in the region \v\ > Ap along with the modified contours Ci 
defined in Appendix C. From the definition of the dressed energy e{v) (|A.7|) and the dressed 



charge Z{v) ( A. 12 ), we have a solution to the above equation. 



Inai(w) = -Pe{v) + ( - 1 ) • 2TTidc + 0(1//?). (4.10) 



Thus we obtain the value 

lnoi(AF) = (^y^-l) ■27ridc + 0{l/f3). (4.11) 
(4) Singlet pair excitations 

As mentioned above, singlet pair excitations are characterized by annihilated 2nc charges. 
Correspondingly, one observes the function 2li(u) has 2nc additional zeros near the Fermi 
points, i.e. ric zeros are located near the left Fermi point and remaining ric zeros are located 
near the right Fermi point. All these additional zeros are distributed in the lower half plane 
and are circumvented by the deformed contour in a manner different from case (3), in counter- 
clockwise manner at both Fermi points (see Appendix C). Simultaneously, ric additional poles 
appear in the upper half plane and are located on the 7r/2-axis. Though these poles do not 
appear explicitly in the NLIE, they lead to a non- vanishing winding number of ai{v), i.e. 
Inai(f) is not periodic (see Fig. 10) 



Im[lnoi(±7r/2)] = T^cVr. (4.12) 

In the low temperature limit, ln2li(f) can be replaced by ln2li(u) ~ for |?;| < Af, ln2li(v) ~ 
Inai(w) — 27rinc for —tt/2 < v < — Af and ln2ti(i') ~ Inai(f) + 27rmc for Af < f < 7r/2 (see 
Fig. QOl). Thus we obtain the integral equation 
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Figure 10: Numerical results for the functions Im In ai{v) (a) and Im ln2li(u) (b) for = 7 = 
1, He = 1.00, /? = 50 and Uc = 1. One observes that the function In ai{v) is a superposition of a 
periodic function and a linearly increasing term. The function ln2li(w) is discontinuous at the 
Fermi points Ap = ±0.748. The discontinuity derives from the additional zeros (circumvented 
by the contour) at the Fermi points. 



In ai{v) 



Ci 



Peo{v) + h * S[\nai]iv)+0{1/P), 



(4.13) 



Ci 



where the symbol * denotes the convolution and the functional S[g] applied to an analytic 
function g{v) is defined stepwise by 



S[g]iv) 



g{v) — 2irinc for —tt/2 < v < — Ap 
for -Ap <v<A-F 

^g{v) + 27rinc for Ap < v < 7r/2 



(4.14) 



with "analytic continuation" out of the interval [— 7r/2,7r/2] by the description S[g]{v + vr) = 
S[g]{v) + Airinc- We directly see that the real part of the solution to ( [4. 131 ) is given by the 
dressed energy ([A. 71 ) 

lnai{v) = -peiv) + if (v), (4.15) 
and the imaginary part / has to satisfy the linear integral equation 

if{v) = h S[if]{v). (4.16) 
Obviously, the solution is unique and odd. Taking the derivative of ( [4.16| ) we find 

nv) = f[>iv)+h%'f'iv), 

f^(v) = ^[h{v - Ap)(/(Ap) + 2nc7r) + ki{v + Ap)(-/(-Ap) + 2nc7r)}. (4.17) 

The symbol Ci denotes the integration contour [— 7r/2, — Ap] |J[Ap, 7r/2]. From the relation 
/(— 7r/2) = — /(7r/2) = ncvr, we find for the odd function f{v) 



2/(Ap) + 2n,7T = - [ f'(v)dv = - [ 
JCi JCi 



z{v)m 



dv. 



(4.18) 
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Substituting ( 4.17 ) and using the definition of tfie dressed cliarge Z{v) ( |A.12 ), we arrive at 



/(^f) = ( 177^ - 1 ) • 27rne. (4.19) 



Z(Af 

From eq. ( [4.15 ) and the above result, Inai(AF) is expressed by 



Inai(AF) = - l) • 2vrmc. (4.20) 

(5) General cases 

Due to hnearity and the results in (l)-(4), one finds 

Inai(AF) = - l) • + (^(^^ " l) ' ^^^^c (4.21) 

4.2 0(1/(3) corrections for the NLIEs 

To evaluate the 0(l//3) corrections for the non-linear integral equations (4.1), we must calcu- 



late the integral term, taking into account its behavior near the Fermi points iLAf. To achieve 
this, we divide the integration term into three parts as 

— / ki(x)ln^2ii(v — x)dx = / ki(v — x)ln^i(x)dx + ricln r- 

vr J_7i J_n cos{v — I'-f) 



1 f cos(v ~\~ ^'y) 

— > — / ki{v — x){lnai{x) + 2TTidc + 2'iTia{x))dx + Ucln -, r 

vr cos{v - Z7) 

+ - [ ki(v-x)ln(l^ ^]dx + -[ kUv - x)ln(l + ai(x))dx, (4.22) 

TTJCi V ^1{X)J J\v\<Af 

with a{x) = —nc,0,nc for x G [— 7r/2, — Af], [— Af, Af], [Af, vr/2]. Note that the additional 
"cos" -term which depends on the definition of the convolution (see Appendix C) appears 
in the above equation. As we have already treated the first linear term in eq. ( [4.22| ), we 
concentrate on the last two non-linear terms. Let us split the integration intervals into two 
parts; [— -7r/2, 0] |J[0, 7r/2]. The integrations over the negative interval read 

-[ ki{v-x)ln(l^ TT]dx + -[ ki{v-x)ln{l + ai{x))dx. (4.23) 

Changing the variable to z = — In ai (v) and using the fact that the derivative of In oi (v) is 
dominated by the one of —j3e{v), i.e. dz/dv ~ — ln'ai(— Af) = /3e'(— Af) = — /^^'(Af), one 
obtains 



ki{v + KY) 



-In ai(-AF) 



7r^e'(AF) yj-lnaiC-Ap) J-oo J 

The integration contours should be taken to circumvent the singularities (see Appendix C). 
The resultant integration is expressed as 

A:i(^; + Af) Ifi^ ai(-AF) + 27ri(4 - nJl^ - 47r2n-V (4.25) 



7r/?e'(AF) V 6 2 
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Adding the results of the positive real axis, we obtain 



— / ki(v — x)ln(l-\ ]—] dx + — [ ki(v — x)lii(l + ai(x))dx 



/3e'(AF) iV6 



4Z\+ ]ki{v - + [-- 4Z\~ ]ki{v + A; 



1 / n 

2 V^f) ^ 2 



±^Z(Af) +n±. 



(4.26) 



Here is nothing but the conformal dimensions defined in ( A.ll| ). Hence we obtain the 
NLIE up to 0(1/(3), 



In oi (v) = -/3e(°) (v) + ki (In ai + 27ridc + 2ma) (v) + Uc In 



Ci /I . ^ • , . ^ • X / ^ , 1 cos(f + Z7) 

cos(f — i^) 

- Mm { (i - - Af) + (i - 44-) k(v + Ap; 

4.3 0(l//9) corrections to the eigenvalues 

Repeating the same argument for the eigenvalues, we obtain 

^({v) + i/j{-v) 



(4.27) 



InA(O) = 2/?ch(a + 1)7 - / 



+ 



'c^ 27rsh(a + 1)7 
27rV'(AF)(l/6 - 2A+ - 2A~) 



(Inai(t') + 2-Kidc)dv 



(4.28) 



/3e'(AF)sh(a + 1)7 

Next we apply the integral identity for bare (oq, 60) aiid dressed functions (a, 6), i.e. if 
a = aQ + k * a, and b = bo + k * b then we have / abo = J aob. Applying this to eq. ( 4.27 ) 
("a = Inai + 27ridc + 27ricr") and the density function ("6 = f?") defined in ( [A.4| ), we obtain 



Z7rsh(a + 1)7 2 



2/3e'(Ap 



^-AA+^ j ki{v - Af) Q{v)dv+ (^^-AA-^ j ki{v + A 



p)Q{v)dv 

(4.29) 



Substituting (4.29) for (4.28) and using ( [A.4 ), we arrive at 



In A(0) = -P{eo + 2{n, - l)ch(a + 1)7 - fin,} - ^'g(AF)(l/6 2^+ 2A ) ^ ^.^^^^ 

pe ( Af ) 

= -p{eo + 2(ne - l)ch(Q + 1)7 - finA + ^(V6 - 2^ - 2Z\ ) ^ ^.^^^^^ ^^ ^^^ 

P^F 

where Eq and wf denote the ground state energy per site (A. 6) and the Fermi velocity ( A.14| ), 
respectively. Using (|4.30| ) and ( 2.29 ), one obtains the low temperature asymptotics of the 
correlation lengths: 



2tt{A+ + A-)' 



(4.31) 



Prom the above calculation and the numerical observation, we can describe the low tem- 
perature behavior of the grand potential ( 2.2^ ) and the correlation lengths ( |2.29| ) by selecting 
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the "quantum numbers" {nc,dc,n'^ ,n^). 
(a) Free energy 

Selecting the quantum numbers (nc, dc, n+, n~) = (0,0,0,0), we obtain the 0(1//?) correc- 
tion for the grand potential per site ( ^.28 ), 



f = {eo + 2K - l)ch(a + 1)7 + f^Ue} - (4.32) 



This asymptotics agrees with CFT with central charge c = 1. 



(b) Density-density and longitudinal spin correlations 

To obtain the low temperature asymptotics of the correlation lengths for the density-density 
correlations {rijUi), the selection rules require to set the quantum numbers to (ric, dc, n~^,n~) = 
(0, ±1,0,0). Thus we obtain 

In addition to the above correlation lengths, one observes a 2kp oscillation term in the corre- 
lation functions 

2/cF = Trree. (4.34) 

In the case of the correlation lengths of the longitudinal spin-spin correlations .^si> we have to 
choose the same quantum numbers as for the density-density correlations. Hence we have 

4i = (4.35) 
(b) Sub-dominant density-density (longitudinal spin) correlations 

One also obtains the sub-dominant correlation lengths for the density-density or longitudi- 
nal spin-spin correlations which are described by the quantum numbers {nc,dc,n~^ ,n~) = 
(0, 0, 1, 0) or (0, 0, 0, 1). The resultant asymptotics of the correlation lengths are given by 

C' = 4T' = ^- (4.36) 
(d) Singlet superconducting pair correlations 

Finally, we consider the low temperature behavior of the correlation lengths ^sp for the singlet 
superconducting pair correlations {cj^i |Cj+ijCj^j^) determined from the quantum numbers 
{nc,dc,n~^ ,n^) = (±1,0,0,0). The result reads 

= (4.37) 

All these analytical calculations of the low temperature properties are consistent with the 
numerical results in the previous section and the predictions from CFT in Appendix A. 
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5 Summary and discussion 



In this paper we have investigated the finite temperature correlations of the strongly correlated 
electron system in one dimension with ?7|j(sl(2|l))-invariance and various types of interaction 
such as Hubbard and correlating hopping terms. As a consequence we observe a competition of 
normal versus "superconducting" correlations. By concrete calculations of the corresponding 
correlation lengths we have found dominant pair correlations for low particle densities n and 
temperatures T. In detail, we have found that singlet pair correlation lengths dominate over 
density-density correlation lengths for low ne and T, but also for high n and T. In the latter 
case, however, the one-particle Green function dominates over all other correlations. 

The computational task was achieved by a first principles approach based on a lattice path 
integral formulation of the Hamiltonian at finite temperature and the subsequent diagonal- 
ization of a suitably defined quantum transfer matrix describing transfer in chain direction. 
The Bethe ansatz equations for the leading and next-leading eigenvalues were reformulated in 
terms of a finite set of non-linear integral equations where the limit of infinite Trotter number 
could be taken analytically. The NLIEs were solved numerically for arbitrary finite temper- 
ature and various particle densities. In the low-temperature limit the analytic investigation 
confirmed the CFT picture and reproduced the dressed energy and dressed charge formalism 
known from the finite size analysis of the T = problem. 

The QTM was based on Z2-grading formulations, which reflect the proper fermionic statis- 
tics. The formulations are applicable directly to other strongly correlated electron systems. 

It is an interesting problem to compare our results for the fermionic system with the ones 
of the corresponding spin system. The eigenvalues of the QTM constructed by the (null 
Grassmann parity) i?-matrix which satisfies the ordinary YBE are written as 

^ ' ' qiiv + ^j) ' ' qi{v + il) Q2iv) 

qi{v-^j) q2{v) qi{v-^7) 

(5.1) 

Due to the extra factor (— )^~" (AT denotes the Trotter number and n is the number of 
electrons), the one-particle Green function for the original system no longer corresponds to 
the one for the spin system. Consequently the characteristic fep-oscillations for the one-particle 
Green function disappear and instead 2A;f oscillations will enter the corresponding correlation 
function in the spin system. 

In our work we have provided the framework for a systematic investigation of the corre- 
lation lengths at finite temperature with an application to a characteristic set of interaction 
parameters. We have not yet performed a comprehensive study of the phase diagrams for a 
wider range of parameters. This remains to be done in a future publication. Also, the effect 
of a non-vanishing magnetic field remains to be explored. We expect that for the isotropic 
limit, i.e. 7 — > 0, the "superconducting regime" in the density-temperature phase diagram 
will shrink to a point. Also the effect of a magnetic field will lead to similar NLIEs where the 
chemical potential terms have to be replaced by /3(/x ± h/2) and 2/3// remains fixed. The 
solutions to these NLIEs are expected to break the symmetries of the Bethe ansatz patterns 
with respect to reflections at the imaginary and/or real axis leading to independent Fermi 
momenta of spin up and spin down electrons. 
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Appendix A Low temperature behavior from CFT 

The ground state properties and the long-distance behavior of some correlation functions have 
already been investigated in |24] with the aid of the root density method and conformal field 
theory. Here we slightly modify the formulation of |24| for direct comparison of our results with 
the prediction from CFT. The ground state properties and low-lying excitations considered 
here are characterized by the BAE 



sin(Aj + (f + 1)70 sin(Aj + f 7^) 



sin(Aj — (f + 1)7*) sin(A 



N,/2 

n 



sin(A/c — \j + i^) 



2 

a 

■J 2 

N 



n 

k=l 



sin(Aj - Afc + ^7) 
sin(Aj - Afc - Z7) ' 



n 



J sin(Afc - Ai + i-j) 



j-J^ sin(Afc - Xj - ij) jj^ sin(Ayt - A/ - i-y) 



(A.l) 



where A'e and denote the total number of electrons and the number of down-spin electrons, 
respectively. We derive directly the above equation from the fact that the ground state energy 
is described by 2-string electron rapidities, and hence shifting the rapidities Xj — > Xj + i'y/2, 
Xjj^i —I- Xj — iy/2 in eq.(6) in |24]. The ground state energy can be written as 

Eo = -2Y^ [cos{k{Xj + i7/2)) + cos{k{Xj - ij/2))] 

i=i 
N,/2r 



2ch(a + 1)7 + 



^(Xj) + 2l;i-Xj) 



where 



k(v) = TT — 2 arctan 



tanv 



th((a + 1)7/2) 



(A.2) 



(A.3) 



By solving the spin rapidities A^ in terms of the charge rapidities Xj, we obtain the density 
function q{v) in the thermodynamic limit, 



Q[V) 



^(w) + ip{-v) 



1 

7rsh(a + 1)7 vr 



ki{v — x)g{x)dx, 



(A.4) 



Ci 



where the Fermi points itAp are determined from the subsidiary condition for the total electron 
density rie 



Q{x)dx = Ug. 



Ci 



Hence the ground state energy per site Eq is expressed in terms of the density function 
.0 = -/^ (2ch(a + l)7 + ^^^^^^)^.(x)dx. 



(A.5) 



(A.6) 
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The dressed energy e{v) 



e{v) = e^^\v) + — I ki{v — x)e{x)dx, 



1 



vr 



(A.7) 



characterizes the low-lying excitations. Here the bare energy function e^^\v) is defined by (as 
in (El)) 

e(0)(?;) = -^^(t,) - ^^(-u) -2//. (A.8) 
Note that the chemical potential is shifted by ^ + 2ch(a + 1)7 /U. The condition 

e(±AF) = 0, (A.9) 
provides another way to define the Fermi points ztAp for given chemical potential fi. 



Density-dGnsity 
Superconducting pair 
2 DensiVdensity (sub) 




Figure 11: Numerical results for the correlations based on CFT (T <^ 1). The density-density 
correlations, singlet superconducting pair correlations and the sub-dominant density-density 
correlations correspond to the quantum numbers {ric, dc,n~^ ,n~) = (0, ±1,0,0), (±1,0,0,0) 
and (0,0,1,0) ((0,0,0.1)), respectively. 

As mentioned in section 4, we observe massless excitations described by quantum numbers 
{nc,dc,n~^ ,n~). Using the concept of CFT, one obtains the correlation functions for primary 
fields as 



(0^±(^)0zi±(O)) 
do: 



exp(— 2idc^F) 



where are the conformal dimensions of the primary operators 



, ^ ± ^Z(Af 
2 V^(Af) 2 ^ ' 



± n^. 



(A.IO) 



(A.ll) 



Here the function Z{Ap) denotes the dressed charge determined from the integral equation 

1 



Z(v) = 2 + - I ki{v - x)Z{x)dx. 
^ Jci 



(A.12) 



For finite temperatures in the scaling limit (0 < T ^ 1), where conformal invariance is 
valid, the correlation functions are obtained from the T = case by the following replace- 
ment § in ([AlO| ), 

VpP TTX 



X 



vr pvp 



(A.13) 
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where vp is the Fermi velocity defined by 

e'{-AF) ^'(Af) . 
VF = = 1--. (A.14) 

TTQiv) Tigyv) 

Thus the long-distance behavior of the correlation functions G{x) is given by 

G{x) ~ cos{2kFdcx) exp ( ^— j , (A.15) 

Fitting this with G{x) ~ exp(— x/,^), one obtains the correlation lengths ^ at T <C 1 



The numerical results of (A. 16) are shown in Fig. 11 for {nc,dc,n^ ,n ) = (0, ±1,0,0) 



(±1,0,0,0) and (0,0,1,0) ((0,0,0,1)). These low temperature asymptotics of the correlation 
lengths are consistent with ( |4.31| ). 

Appendix B Derivation of NLIE for largest eigenvalue 

Here we derive the NLIE for the largest eigenvalue by using the auxiliary functions defined 
by ( |3.lD . The derivation is applicable to the general eigenvalues after a simple modification of 
the integration contours. 

The auxiliary functions are defined by certain combinations of Xj{x) (1 < j < 4): 

Ai(a;)(A3(x) ±A4(x)) ^, , , (Ai(x) ± A2(x))(A2(x) ± A3(x) ± A4(x)) 

A2(x)A(x) A2(x)A(x) 

A2(x) - A2(x) ± A3(x) ± A4(x) 

A3(x)±A4(x) A3(x)±A4(x) 

A2(x) + A3(x) + A4(x) A2(x) ± A3(x) ± A4(x) 

where x = v + ^a'j. Note that the eigenvalues are written in the following form 

A{v) = Xi{v) + X2{v) + Xsiv) + X4{v). (B.2) 
Using the definition of Xj{x), we can write 

A> W + H.) = ^■'-+/^^° + ')) U(.)e-.0 + ^,(„)M^„.,A . (B,3) 

gi(x±§7) t q2{x) J 

One finds that function ( [B.3[) has poles stemming from 52 (x). Due to the BAE ( |2.34| ), the 
poles from qi{x + ^7) must be canceled out. Hence the second factor in (|B.3|) must take the 
following form 

Mx)e''' + Mvf-^^^e^' oc ^^"(^^ + ^^)gi(^ + §7)g2(x)^ 
Q2[X) 0+(x)sh 2 (ix — UN — aj)q2{x) 

iV 

(p+ix) = [sh(zx ± UN ± a7)sh(zx ± un — [a + 1)7)] 2 . (B.4) 
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where the function q2{x) is a certain analytic function which has N + m — n zeros for the 
largest eigenvalue case. Substituting (B.4) for (|B.3|), we obtain 



Ai(x) + A2(x) oc 



N ■ I 

sh~(ix + UN)qi{x + §7(2a + l))q^{x) 



4)j^{x)s\i 2 {ix — UN — 0('y)q2{x) 
Similarly we find A3(x) + A4(x) can be written as 

sh 2 (^ix — UN)qi{x + §7(2a + l))q2{x — ij) 



(B.5) 



A3(a;) + A4(x) oc 



(p^{x)sh. 2 (ix + utv — (a + 1)7)92(2;) 



(/)_(x) = [sh(zj; — UN — a^)sh.{ix — un + (a + 1)7)] 2 . 



(B.6) 



Thanks to the BAE ( p.34D , the eigenvalues ( |B.2| ) are analytic functions (except for the trivial 
poles which derive from the vacuum functions in ( |2.31 ) ) . Thus the eigenvalues ( p. 21 ) must take 
the form 



A(x) oc 



qiiv + ^-f{2a + l))q^{x) 



(B.7) 



where qi{x) is an analytic function having 2A^ — n zeros for the largest eigenvalue case. After 
shifting the parameter x ^ v + ^ctj, we obtain the functional relations to be satisfied by 
certain combinations of the auxiliary functions 



ao(f) oc 



ao(f) oc 



[sh{iv - UN + ^a)sh(w + un - ^{a + 2))]^q2{v + §7(0 - 2)) 
qH^ + io^l)Q2{v + Ha + 2)) ' 

N 

[sh{iv + UN - ^a)sh{iv - un + ^{a + 2))]^q2{v + §7(0 + 2)) 
qiiv + ^jia + lMiv + ^^ia-2)) ' 



2lo('y)2li(u) oc 



(l)+{v + |a7)(/)_(t; + §07) 

qiiv + ^^ia + l))q^iv + ^^a) 



(B. 



2lo(?;)2li(u) oc 



[sh.{iv + un + ^a)sh(it; — un — '^a)\~ q2{v + §7(a + 2)) 

qi{v + ^ia)q2{v + §7Q) 

[sh(ii; + UN + '^a)s\i{iv — un — ^a)]~(72 (^ + §7(o^ + 2)) 



For the largest eigenvalue, we find numerically the following analyticity properties: qi{v + 
|7(a + 1)), q2{v + ^7(0 + 2) (g^ (w + §0:7), (?2 (''^ + §7(0^ ~ 2))) are analytic and non-zero in the 
upper (lower) half plan^, and \2{v + §7) + \^{v + ^7) + A4(u + ^7) is analytic and non-zero 
in the physical strip — ^ < Im t; < ^ . 

Because of these analyticity properties, we can take the logarithmic derivative and perform 
the Fourier transform on both sides of ([B.8D. Explicitly we find 



_^2efc7(«+2) _ E{-l{a + 2) + un) 



UN J 



for A; > 
for A; < 
for A; = 



■* In fact we also observe that qi{v+ fai) = qi{—v+ |7(q! + 1)) and q2{v + §7(0 — 2)) — q2(—v + ^7(0 + 2)) 
are valid in the largest eigenvalue case and some special sub-leading ones. 
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2li 



ai 



2I0 + 2I1 



-9? 



un) for fc > 

i?(-^a + tiAr) for /c < 
for A: = 



qi + 92 

?(^a - uat) - £^(^(0 + 2) - uat) for A; > 

fc7("+i) + + un) + E{-^{a + 2) + ujv) for A: < 

for A: = 



qie 



E{^a + Un) 



qie 



fe7(a+l) 
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fe7(a+2) 



UN) 



2lo + 2ti = < 



1 9^ + 91- 
l9?-^ + 92 



92 



for A; > 
for A; < , 
for A; = 



E{^a + Un) for A; > 
for A; < 
for A; = 



(B.9) 



where we have adopted the notation / for the Fourier transform of the logarithmic derivative 
of a function /(f), 



/ 



d_ 

dv 



lnf{v] 



J2ikv 



dv, 



d_ 

dv 



lnf{v)= f2 



-2ikv 



(B.IO) 



and E{x) := iNe Applying this procedure to the eigenvalue (p.2[), we have 



A 



- E{{a + 1)7 - Un) - E{aj + un) 



A:7(2a+1) 



91 e' 
l9? + 9? 



for A; > 

+ E{-{a + l)j + UN)+E{-aj-UN) for A; < • 

for A; = 



(B.ll) 



From the last two equations in ( p.9| ), we find the Fourier modes qi, and 5^ can be 

expressed in terms of 2to, 2lo and 2li. Substituting them for the first three equations in 
and integrating over v after performing the inverse Fourier transform ( [B.10| ), we arrive at 

Inao(v) = ipNiv) + ko* ln2lo(w) + A;o * ln^^i{v) + Pfi, 
Inao(w) = iPn{-v) + A;o * ln2lo(u) + ko* \d.%i{v) + Pfi, 

Inai(v) = ^Ar(w) + tpNi—v) + ko* ln2lo(u) + ko * ln2lo(w) + ki * ln2li(f) + 2/?//, 



IpNiv) 



N sh.{iv — Un + Za)sh.{iv + un — ^{a + 2)) 

— In . 

2 sh{iv + UN + 3a)sh(it; — un — Z{a + 2)) 



(B.12) 



where the integration contours for the convolutions with InSto, InSlg and InSli are taken by 
the straight line with the imaginary part —6, +5 and {6 is arbitrary but fixed in the range 
< 6 < 7/2), respectively. In the same way, we obtain 



In A""^^(v) = ^Ar(w) + C * ln2to(i;) + C * In^oiv) + (C + C) * ln2li(^;) 

, , N sh{iv + UN + j{a + l))sh{iv - UN - ^{a + 1)) 

^n[v) = — in — — -— . 

2 sh(if — UN + 7(a + l))sh{iv + un — 7(0 + 1)) 

In these NLIE, we can take the Trotter limit — > 00 analytically. 



(B.13) 



lim tpNiv) =^{v), lim ^n{v) = ^(t'). 

N^oo Af— >oo 



(B.14) 
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We have determined the integration constants in (|3.2| ) and in (|3.5|) by taking the 

Hmit I if I oo ^ and comparing the results with 

^o^i^' ^^^2i^' ^-(-'^ + 1)'- 

Appendix C Derivation of NLIE for excited states 

Here we consider the excited states characterizing the correlations which are discussed in the 
main text. As mentioned in section 3, the excited states are described by certain distribution 
patterns of additional zeros and poles which enter the physical strip. To derive the NLIEs, 
we must explore these zeros and poles. By taking the integration contours in a way that they 
circumvent the zeros in the physical strip, we can derive the NLIEs directly by utilizing the 
above results for the largest eigenvalue. 

Taking the integration contours of the convolutions in ( [B.12| ) as in Fig. |l^ (and subsequent 
figures), we have the same algebraic structure of NLIEs as in ( p.l2| ) without additive terms. 



This is true at least in the case of the intermediate (differentiated) version of (B.12) where 
In a. In 21, ipN are replaced by (In a)', (In 21)', ip'j^. 

When integrating the NLIEs with respect to the argument we want to make sure that the 
convolutions are understood in the manner k * (ln2l)'(t') := /c(x)(ln2t)'(f — x)dx with a 
i;-independent path C. This can always be achieved with C starting (ending) at —tt/2 (+7r/2) 
and suitable deformations in order to avoid the singularities. The integration with respect to v 
and up to an integration constant yields a set of NLIEs which has the same algebraic structure 
as in the case of the largest eigenvalue where now terms k * ln2l(?;) = k{x)ln^{v — x)dx 
appear with deformed contours. Alternatively, we can straighten the contours and due to 
Cauchy's theorem obtain additional terms from the singularities of the auxiliary functions 
(these are the "sin"-terms below). 

On top of these terms we obtain a second set of additional terms. Note that the functions 
hi^oiv) and \d.%i{v) are not periodic: ln2lo(7r/2) - ln2lo(-7r/2) = 2m and ln2li(7r/2) - 
ln2li(— 7r/2) = 27ri for the one-particle Green functions and ln2lo(vr/2) — ln2lo(— 7r/2) = 47ri 
and ln2li(7r/2) — ln2li(— -7r/2) = 47ri for the triplet superconducting pair correlations. Further 
note that A;*ln 21(11) = k[x) \x].'^[v — x)dx = k{v — x) ln2l(3;)(ix which in the general case 
(non-periodic In 21) is different from J_j^ k{v — x) ln2t(x)(ix. For many numerical applications 
we actually like to have convolutions of the last type where the auxiliary functions are evaluated 
on fixed (straight) contours. The difference of the two versions of convolution integrals can be 
worked out and gives rise to the "cos" -terms below. 

By use of the analytic continuation of the NLIE (see Appendix D), one determines the 
integration constants in a manner similar to the largest eigenvalue case. 

We want to categorize the additional zeros and poles of the auxiliary functions in the 
following four cases of excitations: (1) the one-particle Green functions and the triplet su- 
perconducting pair correlations, (2) the density-density (longitudinal) correlations, (3) the 
transversal spin-spin correlations and (4) the singlet superconducting pair correlations. 

(1) One-particle Green function and triplet superconducting pair correlations 

First we consider the NLIE for the one-particle Green functions (cj^^Cj^o-)- In this case, the 
eigenvalues ( |2.31| ) belong to the sector n = N — 1 and m = N/2 — 1. From the numerical 
calculations with finite Trotter number A^, we observe additional zeros, namely a zero of 

^ Note that we take this hmit by making use of analytic continuation of the NLIEs (see Appendix D). 



28 



(a) 



Im V 



(b) 



Imv 





Co 


Co 




Ci 


^Rev _| 


Ci 




Co 








o: 6>i 


o: 6*1 □: -e* 



9o 



Figure 12: Integration contours for (a) the one-particle Green function and (b) the triplet 
superconducting pair correlations. Here Co, Co and Ci denote the contours for 2lo, 2to and 
2li, respectively. 



the auxiliary function 2lo(w) and a zero 9i of located in the lower half plane. The zeros 

^0 and 01 are derived from 

i \ ./„ I \ ./„ I \ „ h / " * 



A2(^o+27al +A3Uo+27«j +A4Uo+27aj =0, gM0i + -7a)=O. (C.l) 
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Figure 13: Integration contours of the NLIE for (a) the dominant and (b) sub-dominant 
contributions to the density-density correlations. 

In the case for the triplet superconducting pair correlations {cj_^_i ^cj |Ci+i^|Cj^|), the corre- 



sponding eigenvalues ( 2.31 ) describing these correlations belong to the sector n = N — 2 and 
m = N/2 — 2. Numerically we observe additional parameters: zeros Oq, —9q for the function 
2lo(t') and zeros 9i, —91 function 2ti(w) appear in the lower half plane (the symbol * denotes 
complex conjugation). As for the one-particle Green function, the zeros ^o (~^o) ^^'^ (^i) 
satisfy the first and the second equation in ( p.l| ), respectively. 

(2) Density-density correlations 
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Figure 14: Integration contours of the NLIE for the transversal spin-spin correlations. 

Let us consider the NLIEs for the density-density correlations {rijUi). The eigenvalues ( |2.31| ) 
belong to the sector n = N and ra = N/2 which is the same sector as for the largest eigenvalue. 
These correlations are described by the distribution patterns of a zero 6*0 (— ^'o) and a pole 
di of the auxiliary function 2lo(w) (2to(f )) appearing in the physical strip. In this case, 

these zeros and poles are derived from 

<i (Oo + ^7") = 0, q2 (^-Oo + ^7a) = 0, 

Qi (Oi + '-ja^ = 0, qi (^-Oi + '-j{a + 1)^ = 0. (C.2) 

In addition, we can consider the sub-dominant terms described by simple particle-hole exci- 
tations at each Fermi point (see section 4). These sub-dominant terms are determined from 
zeros ^0 (^o) 2to(f) {^o{v)) and 9i, 9^ for 2ti(w), respectively. These zeros and poles satisfy 



eq. (C^) (replace —9i and — by —9i — > 91 and —9o — > 9q, respectively). Taking the inte- 



gration contours as in Fig. 13, we arrive at NLIEs which have the same form as in ( p.l2| ). For 
both cases, we take the convolutions in the ordinary way, i.e. as for the largest eigenvalue. 

(3) Transversal spin-spin correlations 

Let us consider the eigenvalues which describe the transversal spin-spin correlations {(t^(t^)- 
The eigenvalues ( p. 31 ) belong to the sector n = N and m = N/2 — 1. Making use of numerical 



calculations, we observe additional zeros (~^o) and poles 9i {—9i) of the auxiliary function 
2to(u) (2lo(f)) appear in the lower (upper) half plane. As in case (1), these zeros and poles 
get close to but never cross the real axis in the low temperature limit. These additional zeros 
(poles) are derived from 

Aa (^±9o + ^7a^ + A3 (^±^0 + ^7"^ + A4 (^±9o + ^70^ = 0, 

Qi (^1 + ^7«) = 0, qi (^-9i + ^7(a + 1)) = 0. (C.3) 

We derive the NLIEs by taking the contours as in Fig. U. 

(4) Singlet superconducting pair correlations 
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Figure 15: Integration contours of the NLIE for the singlet superconducting pair correlations 
for T > Tc (a) and for T < (b). 



We consider the singlet superconducting pair correlations (c^_)_]^|''^ 



|Cj+i^|Cj_l). In this case, 

the eigenvalues written as (2.31) belong to the sector n = N — 2 and m = N/2 — 1. These 
correlations are described by the distribution patterns of zeros Go-, do and poles 6i, —Of of 
the auxiliary function 2to(u). Numerically we observe that the zero Oq is located on the ib7r/2- 
axis (Ref = ±tt/2) and moves to the real axis with decreasing temperature. At sufficiently 
high temperatures, the zeros and 9q are related by = close to the boundary 

(ib7r/2-axis) with decreasing temperature. At a certain temperature Tc, the zeros 9o, Oq and 
^0 satisfy Oq = 9o = Oq- Below Tc, the relation = — longer valid, and these zeros are 
located on the 7r/2-axis satisfying the relation Im^o < Im^o < Im^o- and move upwards 
(into +i direction along the 7r/2-axis), opposite to this moves into —i direction. In the low 
temperature limit, comes near to the real axis but never crosses it. On the other hand 
enters the upper half plane and finally leaves the physical strip. Oq also leaves the physical 
strip in the lower half plane. Consequently at low temperature, only 6i and —9\ (zeros of 
%i{v)) are circumvented by the integration contour of 2ti(f) and these zeros characterize the 
massless pair-excitations (see section 4). 

Due to this behavior, we must take different integration contours for the different regimes 
T > Tc and T <Tc. These zeros are derived from 



A2 ( (^o + ) + A3 ( ^0 + j + A4 ( ^0 + ^7" 



0, 
0, 



q2{Oo) =0, q'l{ei + -7a j =0, qi( -01 + ^7(q + 1) 



0. 



(C.4) 



We derive the NLIEs by taking the integration contours as in Fig. 15 a for the case T > Tc and 



as in Fig. |l5| b for the case T < Tc 



Additional terms in the NLIEs 



As mentioned above, we obtain the same structure of NLIEs as in the case of the largest 
eigenvalue by taking the integration contours of the convolutions in a way that they circum- 
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vent the zeros. To solve the NLIEs numerically, we are forced to use straight integration 
contours. Due to Cauchy's theorem in this formulation we obtain certain additional terms in 
the NLIEs. In the remainder we list the additional contributions to the driving terms in (|3^ 
and (|1 



(i) One-particle Green function 

First we consider the one-particle Green function. The corresponding NLIEs are characterized 
by the zeros and Oi. In the low temperature limit, these parameters come near to the real 
axis but never enter the upper half plane. From (|C.l| ), we find that these parameters satisfy 
the subsidiary conditions 

2lo(eo) = 2lo(^o) = 0, 2lo(0i)=oo, 

2li(^i) = 0, 2li(0o) = oo. (C.5) 

Performing the Fourier transform and employing Cauchy's theorem, we determine the 
NLIE ( ^.71) with the following additional terms 

for Im^n < 



_(0), , I cos{»;+i<5) 
V'O (^) = \ ( n ■ , 



for -5 < Im6'n < 



— lor —o<^ nnun 



cos{v+i&) 
(1), , J cos(v— 17+4(5) 

, sin{v—6\—i^) 



for Im^i < —6 
for -5 < Im 6*1 < 



_ I . _(0), X , . su).{v - Oi + i-f) I \- / \ /n «^ 

= ip'Q'(v)ipl'(v), ipQ{v) = — — — , ipi{v) = ipQ{v)ipQ{v). (C.6) 

sml^u — u\) 



The corresponding eigenvalues h.{v) are written as ( |3.8| ) with the term x(^) (split into two 
terms xi'^) = X^^\'^)x^^\'v)) ^ 



cos(v+ ^■ya+iS) 

X^'Hv) 



cos(i)— ^7(o+2)+i(5) 

sin(v— 60 + 570) 
sm(t)-eo-|7(a+2)) 



(1), X sm{v - 6*1 - §a7) 

'(V) = : X < 

sin(t; - 01 + §7(a + 2)) ^ 



for Im^o < ~^ 
for -5 < Im6'o < 



cos(f + ^7Q+ii5) 



for Im^i < —5 

(C.7) 



sin(D— 6*1 — ^7(0+2)) „ [■ T /I 

— . , \ — for -5 < Im6'i < 

sm(i;-6»i + ^a7) 



In fact the NLIE (3.7) with the additional terms ( p.6| ) have two solutions. These two solutions 
are related to each other by taking the mirror image with respect to the imaginary axis. Cor- 
respondingly, the eigenvalues A(0) are doubly degenerate in magnitude and they are complex 
conjugate to each other. The oscillatory behavior of the one-particle Green function (referred 
to as /cp-oscillation at zero temperature) derives from this degeneracy. The correlation length 
for the one-particle Green function is given through formula ( 2.2S| ). 

(ii) Triplet superconducting pair correlations 

We determine the additional terms for the triplet superconducting pair correlations. As men- 
tioned above, these correlations are characterized by the additional zeros and 9i. In the low 
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temperature limit, these zeros come near the real axis but never enter the upper half plane. 
From (C.l), one finds that these zeros satisfy 



2lo(^o) = M-Oo) = 2to(^o) 
2ti(ei) =2li(-0i*) = 0, 



0, 2lo(^i)= 2lo(-0i*) = 
2li(^o) = 2ti(-^o*) 



(C.8) 



Performing the Fourier transform and employing Cauchy's theorem, we determine the 
NLIE (p.7[) with the following additional terms 



cos^ {v—i-/+iS) 
cos^Jv+iS) 

sm{v—8o —i'y) sin(?j+6'Q —if) 
sin{v—9o) sin(D+6'Q) 

cos^ (v+iS) 
cos-^ (v—i-/+iS) 

sin{v—9i) sin(D+6'j ) 
sin(i;— Si —i'y) sin{v+9'^ —if) 



(pi{v) = ipo{v)lfQ{v). 
The corresponding eigenvalues A(u) are written as ( |3.8[ ) with the term x{v 

cos'^ {v+ ^fa+iS) 



for Im ^0 < —S 
for -6 < ImOo < 
for Im^i < —5 

for -6 < Im^i < ' 

sin(?; — ^1 + i'y) sm(v + 6i + i'y) 
sm{v — 9i) sm{v + 91) ' 



X^'\v) 
X^'\v) 



cos^{v-^f{a+2)+i5) 

sin{v—do + ^fa) sm{v+dg + ^'ya) 
^ sin(i;-6lo-|7("+2)) sm(v+e*-^'y{a+2)) 



for Im^o < ~^ 
for -S < lm9o < 



sin(t) — 6*1 — |a7) sin(u + 91 — |a7) 
sin(7; -9i + ^-f{a + 2)) sin(w + 0* + ^-f{a + 2)) 



' cos'^{v-^f(a+2)+i5) 



X < 



cos^ (v+^-fa+iS) 



sin{-u- 



-|7(a+2)) sm(t)+e*-i7(a+2)) 



sin(ii— 9i + |q!7) sin(D+Si + 507) 



for Im^i < —5 
for -(5 < Im6'i < 



(C.9) 



(C.IO) 



(C.ll) 



(iii) Density-density (longitudinal spin-spin) correlations 

Next we determine the NLIE for the density-density correlations. As the critical exponents 
of the density-density correlation functions are identical to those of the longitudinal spin-spin 
correlations ((7|ct?), the resultant NLIEs are also identical. From ( p. 21) , we find these zeros 
satisfy the subsidiary conditions ( p.l2| ) . 

%{9o)=M-9o) = 0, 2lo(0i)=2to(-ei) = oo, 2ti(±0i) = 0. (C.12) 

Correspondingly, the additional terms in the NLIE ( |3.7D read 



smiv-eo)cosiv-if+iS) _^ 
_(0) / N I sm{v—0Q—if)cos{v+i5) ^ 



1 for -6 < lm9o 

cos(v+i5) 



, , sin(t; + ^1 - ij) cos{v-if+i5) for Im 9^ < -5 
sm{v + 9i) for-5<Im0i<O 

V sm(v—Oi—tf) ^ 
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The additional term x(^) i^i (|3.8D is determined from 



(C.13) 



X^'Hv) 



AO) I 



' sin(ti—6o — 57(0+2)) sin(ti+6o + f 7(0+2)) cos(v+ ^ya+iS) cos{v— ^'ya—iS) 
sin(D— 50 + 570) sin(i)+9o — 570) cos{v— ^■y{a+2)+iS) cos(i)+^7(o+2)— i(5) 



1 



for Im^o 
for -6 < 



< -6 



sm{v — 61 — ia^y) sm{v + + f 07) 



X < 



sin(?; - 6*1 + §7(0 + 2)) sin(f + 6*1 - 

' cos(i;-|7(a+2)+i<5) cos{v+^-y(a+2)-iS) 
cos{v+ ^ja+iS) cos{v— ^-ya—iS) 

sin(t)-ei-i7(a+2))sm(^+6>i + i7(a+2)) 
sin(i'— 61 + ^07) sin(i)+9i — ^07) 



§7(a + 2)) 
for Im^i < —6 

for -5 < Im^i < 



(C.14) 



In the low temperature limit, 9i gets close to the real axis but does not cross it. In contrast, 
^0 crosses the real axis and enters the upper half plane. Hence only one additional zero 
61 appears explicitly in the above NLIE, which is characteristic of massless excitations (see 
section 4). The solutions to the above NLIE have the symmetry: Oo(w) = Oo(— f) and ai{v) 
is symmetric with respect to the imaginary axis. The NLIEs ( |C.13D have two solutions which 
are connected with each other by taking the mirror image with respect to the imaginary axis. 
Hence the corresponding eigenvalues A(0) are degenerate in magnitude; they are complex 
conjugate to each other. Due to this degeneracy, the correlations have an oscillation factor. At 
zero temperature, this oscillatory behavior is the characteristic 2A;F-oscillation in the density- 
density correlations. 

The sub-dominant terms discussed above are characterized by another distribution pattern 
of zeros and poles. Consequently, a zero Oq (Oq) and a pole 9i of the functions ao{v) (ao{v)) 
appear in the physical strip. From ( C.2D (replace —9i and —60 by —9i — > 9'^ and —9o — > 9q, 
respectively), we find that these zeros and poles satisfy the subsidiary conditions 



0(C0J 



21, 



01^0 J 



0, 2t, 



01^1 J 



21, 



0(^1. 



00, 2ti(0i) =2li(0^) =0. 



(C.15) 



In the same way as before one obtains the NLIEs characterized by the following additional 
terms 



' sin{v—9o) cos{v—i'y+i5) 



for Im^n < —S 



sin(ti— 60— *7) cos(D+j(5) '^'^'^ ^^'^'^ ^ 

1 for —5 < Im( 



Sin(t' — 91 — i-y) J cos{v-i-y+i5) 



sin(f 



sin{v—6i) 
sin(D— Si —17) 



for for —5 < Im^i < 



foiv) = (Po\v)ipl^\v), ipo{v) = V*o{v), (pi{v) = ipo{v)ipQ{v). 



(C.16) 



The term x(^) the eigenvalues ( ^.81) is written as 

xiv)=X^'Hv)x^'\v), 

X^'Hv): 



sin(i;- 



■|7(a+2)) sin(D-6'g + |7(a+2)) cos{v+ ^ya+iS) cos{v- ^-ya-iS) 



sin(v— ^0 + 570) sin(i)— 0Q — ^70) cos{v— ^'y(a+2)+i5) cos{v+ ^■y{a+2)—i5) 



1 



for Im^o < 



for —6 < Im( 
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sm{v — 01 — |a7) sm{v — 91 + §07) 
sm{v - 6*1 + §7(a + 2)) sm(t; - 6^* - §7(0 + 2)) 

' cos(ii-|7(a+2)+j<5) cos(-u+57(a+2)-i<5) 



cos{v+ ^^a+iS) cos{v— ^-ya—iS) 
s m(ii-6>i-i7(a+2))sm(i>~q + f7(a+2)) 
sin(ii— 6i + iQ!7) sin(?j— — ^07) 



for Im^i < —6 
for -6 <lm9i < 



(C.17) 



At low temperatures, the parameter 9i (zero of the function )) moves toward the real axis 
but never crosses it, whereas Oq (zero of the function 2to(f)) crosses the real axis and enters 
the upper half plane. In this limit, the NLIEs depend only on the two zeros (of 2ti(w)) 9i and 
91, which characterize the particle- hole excitations (see section 4). 

Note that the solutions to the above NLIEs satisfy a.Q{v) = ao(f) and ai{v) is symmetric 
with respect to the real axis but not to the imaginary axis. Therefore ai{v) is real for f G M. 
The three NLIE are also reducible as in the dominant case. The NLIE with ( |Cl6| ) admits two 
solutions. In contrast to the above dominant case, they are related to each other by taking the 
mirror image with respect to the real axis. The corresponding eigenvalues A(0) are degenerate 
and are real. Hence one finds that the corresponding correlations do not have oscillatory terms. 

(iv) Transversal spin-spin correlations 

Let us determine the additional terms of the NLIEs describing the transversal spin-spin cor- 
relations. As mentioned above, these correlations are characterized by the zeros and poles 



satisfying eq. (C.3). Hence we find that these zeros (poles) satisfy the following subsidiary 
conditions. 



2lo(±^o) 
2li(±ei) 



2lo(±eo) 
0, 



0, 2lo(0i) = 2lo(-0i) = 00, 
2ti(±6lo) = 00. 

Consequently, the additional terms are written as follows: 

for Im^o < ~^ 
for -6 <lm9o < 

for Im^i < —5 

sm{v—6i) 



(C.18) 



cos(d— i7+i(5) 
cos(i'+i(5) 




cos(t)+i(5) 
cos(d— i7+i<5) 



sin(t> + 9i) 



sm{v~9i —27) 

^o(^) = ^o\^)^o\^)^ Mv) = ^o{-v), 



for -5 <lm9i < 
fiiv) = fo{v)ipQ{v). 



(C.19) 



The term xi'^) iii the eigenvalue (|3.8|) is written as 



Xiv) 



X^'\v) 



X'-'{v)x^''[v) 

cos(v+ ^■ya+iS) cos{v— ^•ya—iS) 
cos(v— |7(a+2)+j(S) cos(»;+|7{o+2)— i(5) 

sin(u— (?o + f 70) sin{v+9o — ^^a) 

sm{v — 9i — §07) sin(t> + 9i + |a7) 
sm{v -9i + i7(a + 2)) sin(?; + 9i - ^-f{a + 2)) 

' cos{v- ^'y{a+2)+i5) cos{v+^'y(a+2)-iS) 



for Im ^0 < 
for -5 < Im^o < 



cos{v+ ^ja+iS) cos(v— ^^a—iS) 
s m(t;-ei-|7(a+2))sin(^+6>i + l7(a+2)) 
sin(i;— 6i + ia7) sin{v+Oi — ^a-y) 



for Im^i < —6 
for -5 <lm9i < 



(C.20) 
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As in the largest eigenvalue case, the solutions Oo(w) and Oo(w) have the symmetry ao(f ) = 
Oo(— v). Therefore one reduces the above set of three NLIEs to the set of only two NLIEs. 
In addition, one finds that the NLIEs ( C.19| ) admit two solutions which are related to each 
other by taking the mirror image with respect to the imaginary axis. Thus the corresponding 
eigenvalues A(0) are degenerate in magnitude and are complex conjugate. Hence, oscillatory 
behavior of the correlation functions is observed, referred to as 2A:F-oscillation for zero tem- 
perature. 



(v) Singlet superconducting pair correlations 

Finally we consider the NLIE for the singlet superconducting pair correlations characteriz- 
ing the above mentioned additional zeros (poles). From ( C.4 ), one sees that these zeros and 
poles satisfy the subsidiary conditions 

2lo(^o) = 2to(^o) = 2lo(^o) = MOo) = M9o) = 0, 
2lo(^i) = 2to(-0i*) = 2to(^i) = M-ei) = %i9o + i-f) = oo, 
2ti(^i) = M-oi) = 0, M0o) = MOo) = oo. 

The additional terms in the NLIE (| 



sm{v—6o —i'f) cos{v+i5) 
sin(t)— 6o) cos (u— 17+1(5) 

cos(i)+«(5) cos(t)— 17) 
cos(i)— «7+i(5) cos{v) 

sm{v—8o) cos{v—i-f+iS) 

sin{v—6o) 
sin{v—9o—i'y) 



) read explicitly, 

for Im^?o < —S 
for -6 < Im6lo < , 
for < Im 6*0 
for Im^o < —S 
for -6 <lm9Q < 
for Im^o < ~S 



E2±±i^ {oi-5<lmeo<0 



cos(i)— 47+1(5) 
1 

sin{v—9i) sin{v+9^) cos'^ (v—i-f+iS) 
sin(ti— 61 —47) sin(u+6J' —i^) cos^ (v+iS) 



for Im 6*0 < 



for Im^i < —6 
for -5 <lmdi < 



for <lm6o < 6 



1 

sin(t)— 0o+*7) cos(d) 
sin(t!— 60) cos(tJ+i7) 

cos(^v) cos{v+^^-^5) fo^ 5 < Im 
^ cos(i)+«7) cos(v— 10) ^ 

■ sin(.-eo-n) forIm^o + 7<'^ 

sin(D— 60) 

for5<Im^o + 7' 

. cos(d+27— 1(5) " ' 

sin(T; — 61 + i'y) sm{v + 91 + i^) 



(C.21) 



(0)( 



sin{v — 9i) sin{v + 91) 



(C.22) 
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The additional term x(^) i^i 



is determined from 



X^'Hv) 



x'^'Hv) 



sin{i)— 60 + ^70) cos{v— §7(a+2)+«(5) 
sin{?j— 60 — §7(a+2)) cos(v+^"fa+i5) 

sin(i'— So — §70) cos (i)— §7(0+2) +i(5) cos(ti+§7o) cos(t;+ §7(0+2)) 
sin(D— So + §7(0+2)) cos(v+§70+i(5) cos(ti— §7(0+2)) cos{v— §70) 
cos(ii— §70— j5) cos(i;— §7(0+2) +i6) cos(v+§7o) cos(r!+§7(o+2)) 
cos(ri+§7(o+2)— i(5) cos(ti+§70+j<5) cos(i;— §7(0+2)) cos(i'— §70) 

1 for Im^o < 

sin(i;— So+§7o) cos(i;— §7(o+2)+j<5) 



for Im^o < 
for -6 < ImOo < 
for < Im 00 < ^ 
for 6 < Im^o 



^ sin(D—So — §7(0+2)) cos(f+§70+i<5) 
(1 

sin(i;-So-§7(o+2)) cos(i;+§7(o+2)-i<5) 

sin(i;— 6o + §7o) cos(v— §70— i5) 
cos(v-§7(o+2)+i5) cos(-!;+§7(o+2)-i(5) 



cos(i)+§70+i<5) cos(t;— §70— i<5) 

sin(t> — 61 — §07) sin(z; + 9\ — §07) 
sin(z; -61 + ^j{a + 2)) sin(t' + 0J + §7(0 + 2)) 



for -5 < Im6'o < 

for < Im ^0 + 7 < <5 
for Im ^0 < —S, (5 < Im 00 + 7 
for -5 < Im^o < 



X < 



sm(t;-Si- §7(0+2)) sm(i;+S*- §7(0+2)) cos'^ (v+^-fa+iS) 
sin(D— Si + §07) sin(-t;+e* + §07) cos2(v— §7(o+2)+j<5) 



for Im^i < —5 
for -6 < ImOi < 



(C.23) 



Appendix D Analytic continuation for NLIE 

The above derived NLIEs are weh-defined only in the region —7 + 5 < Imv < for Inao(T^), 
< Imv < "f — 5 for Inao(v) and —5 < Imv < 5 for Inai(v). To determine the integration 
constants (which should be determined from the limit i\v\ 00) and the additional zeros 
for the excited states, we must extend the region by using the analytic continuation. We 
can achieve this by applying Cauchy's theorem to the convolutions. As a consequence, the 
additional terms appear due to the residues of the kernels. The resultant equations can be 
written explicitly as follows. 

In 00(7;) = ip{v) + ko * ln2to(t;) + ko * ln2li(T;) + PfJ,- (foiv), 



lnao(f ) = ip{-v) + ko * ln2lo(u) + ko * ln2ti(u) + f3^i- Tp'o, 

Inai(f) = ip{v) + Tp{-v) + ko ln2lo(u) + ko ln2lo(?;) + ki ln2li(t;) + 2/3/x - ipl{v), 

(D.l) 
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where 



ln2lo(u + ij) + ln2li(v + i^) for Imv < —7 



for —7 < Im V < —7 + S 
for —7 + (5 < Im V <0 
for < Imv < 6 
for (5 < Imf 



ln2lo('u + i'j) 
ifliv) = <j 

ln2ti(w) 

_ln2lo(f) +lnS2li(t;) 

ln2lo(u) + ln2li(u) for Imv < —5 

In 211(f) for -6 <lmv < 

^o(^) = ^0 forO<Imv<7-(5, 
ln2lo(t' — i'j) for 7 — (5 < Im?; < 7 

\n^2Lo{v — ij) + ln^i{v — ij) for 7 < Imv 

ln2lo(v) + ln2lo(v + ij) + ln2li(f + ij) for Imv < —7 

ln2lo(v) + ln2lo('y + ^7) for —7 < Imv < —7 + (5 

ln2lo(v) for —7 + (5 < Imv < — (5 

</?5(v) = .{ for -(5 < Imv < (5 . (D.2) 

ln2lo(v) for (5 < Imv < 7 — (5 

ln2lo(v) + ln2to(v — 27) for 7 — 5 < Imv < 7 

ln2lo(v) + ln2lo('y — ij) + In2li(v — ij) for 7 < Imv 

Similarly, the NLIE for the eigenvalue can be written as follows. 

InA(v) = *(v) + C * ln2lo(v) + C *ln2to(i^) + (C + C) * ln2li(v) + Xoiv) + Xoiv) + xUv) , 

(D.3) 

where 



roiv) 



ln2lo(v + 1 70;) 




for Imv < —^a — 5 



for 



5 < Imv < ^{a + 2) - S , 



_ln2to(v - §7(a + 2)) for ^{a + 2)-5< Imv, 
ln2lo(v + i7(a + 2)) for Im v < -^(a + 2) + 5 



for 



r(a + 2) + (5 < Imv < + S , 



for + 5 < Imv, 



Xliv) 



-i(a + 2) 





_ln2lo(v - §70;) .w. 2 

ln2li(v + §70) + ln2li(v + 17(0 + 2)) for Imv < 

ln2ti(v + |7a) for 

for 

ln2ti(v — §70) for ^ 

^ln2li(v - ^ja) + ln2li(v - §7(a + 2)) for ^(a + 2) < Imv. 



-^{a + 2) < Imv < - 
— ^a < Imv < 
for < Imv < ^(a + 2) 



(D.4) 
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